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We propose that a pseudospin ferromagnetic (i.e. inter-wire coherent) state can exist in a system 
of two parallel wires of finite width in the presence of a perpendicular magnetic field. This novel 
quantum many-body state appears when the inter wire distance decreases below a certain critical 
value which depends on the magnetic field. We determine the phase boundary of the ferromagnetic 
phase by analyzing the softening of the spin-mode velocity using the bosonization approach. We 
also discuss signatures of this state in tunneling and Coulomb drag experiments. 
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Ferromagnetism (FM) in low dimensional itinerant 
electronic systems is one of the most interesting subjects 
in condensed matter physics. As early as in the '60s 
Lieb and Mattiai (LM) has proved that a ferromagnetic 
state cannot exist in one-dimensional (ID) system if the 
electron-electron interaction is spin/velocity- independent 
and symmetric with respect to the interchange of elec- 
tron coordinates. Therefore, possible candidates for ID 
FM must involve some nontrivial modification in the 
band structure and interaction to avoid the restrictions 
of LM's theorem. Most of the examples proposed in the 
literature^ rely on some highly degenerate flat bands (or 
at least systems with the divergent density of states) and 
can be understood as a generahzation of Hund's rula^. 
The only exception appears to be a model of finite range 
hopping with a negative tunneling energy^. 

From the experimental point of view, however, phys- 
ical realization of the ID FM in thermodynamical limit 
is still absent to the best of our knowledge. In two- 
dimensions (2D), some of the most intriguing ferromag- 
netic systems are the quantum Hall (QH) bilayers at the 
total filling factor one. In these systems the flat band 
structure is provided by the magnetic field (Landau lev- 
els) and clear experimental evidence of the 2D pseudospin 
ferromagnetism (PSFM, with the pseudospin being the 
layer index) has been observed in the tunnelingi' and drag 
experiments^ several years after theoretical proposals^. 

In this paper we propose a realistic one- dimensional 
system which should exhibit a pseudospin ferromagnetic 
order. The system consists of two finite-width quan- 
tum wires with a magnetic field applied perpendicular 
to the wire surface, see Fig. ^a). In the presence of the 
magnetic field, single-electron states are modified, which 
leads to a strong effective mass enhancement and modi- 
fication of the effective Coulomb interaction. These two 
effects can lead to a softening of the spin mode veloc- 
ity when the inter-wire distance becomes smaller than 
a critical value dc- The system then becomes an easy- 
plane PSFM state due to the appearance of interwire 
coherence (IWC), which should manifest itself in the ap- 
pearance of the resonant peak in the tunneling conduc- 
tance at small bias voltages. We also calculate the drag 
resistance of such ID PSFM state within the mean-field 
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FIG. 1: (a) Schematic double wire system considered in this 
paper, (b) and (c) are single particle energy -EjJ(fc) of each 
wire for B — and B cases respectively. 



approximation and demonstrate that the drag resistance 
first increases and scales with the longitudinal size as 
the magnetic field is increased (or the inter-wire distance 
is decreased) toward the phase transition boundary and 
then becomes dramatically reduced (i.e. not scaled with 
the size) when entering the PSFM state. The proposed 
ID PSFM transition should be experimentally accessible 
by the present or near future semiconductor technology. 

The double wire system we consider is aligned in the y 
direction, Fig.^a), and centered at a: = and z = ±d/2. 
Electrons are confined by a parabolic potential, imwgcc^, 
in the x direction and their motion in z direction is as- 
sumed to be totally quenched. Using the Landau gauge, 
the single particle Hamiltonian of momentum k in each 
wire can be derived to be 
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where m* = m{uj'^ 
tron mass, oj = \/uj^ 



ll!q)/ll>q is the renormalized elec- 
'2 is the Landau level split- 



/pfc is the guiding center coordi- 
ujc/'m{Lu^ -\- lOq) being the magnetic 



ting, and xq 
nate with Iq 

length. ujc — eB/mc is the bare cyclotron fre- 
quency. The wave-functions and energy spectrum of 
Eq. are easily found from the analogy with the 
standard QH system^: 'iljn,k,sii^) = L^^^'^e^^^ ipn{x + 
xo)^/5{z- sd/2) and El{k) = (n + \)Cj + where 
n is the Landau level index and s — ±i is the 
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FIG. 2: Calculated critical interlayer distance, dc, as a func- 
tion of magnetic field (cx uic)- Electron density in an indi- 
vidual wire. He, is 0.6, 0.7, • • • , 1.0 x 10^ cm^^ from top to 
bottom. Here Ao — 500 A and luo ~ 0.05 meV. 



pseudospin index for the upper/lower wire, ipnix) = 
(7ri/22"n!fo)-i/2g-xV2/^^^(2,/[p) ig ^l^^ eigenfunc- 

tion of a parabolic potential with Iq = -^/l/mo). 
Throughout this paper we concentrate on the strong 
magnetic field (or low electron density) regime so that 
only the lowest energy level {n — 0) is occupied. One 
can see that the magnetic field not only modifies the band 
splitting, but also increases the effective mass in the lon- 
gitudinal (y) direction, leading to a flatband structure 
with high density of states similar to the Landau level 
degeneracy in 2D system, see Fig. ^b). 

The interaction Hamiltonian can be derived to be^: 
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where 

Q.±_ = LW is the wire area. 



electron field operators, 
and T4,s'(q±, fe) = 

«92<i _ 2)j e~«9a;(fei-fc2)'o 

ID interaction with ^(q) be- 
interaction. The form-function, 

is obtained by 



is the effective 
ing the Coulomb 

A(q_L) = exp 

integrating the electron spatial wave function**. Due 
to the presence of magnetic field, the effective ID 
interaction, V^s.s' (qj- 1 ^i; ^2), is not equivalent to any 
spin-independent (or velocity-independent) symmetric 
potential. Thus in our system, the ferromagnetic state 
is not inhibited by the LM's theorem. 

Starting from Eqs. Q-©, one can use the standard 
bosonization approach to describe the low energy physics 
near the Fermi points. After neglecting the irrelevant 
(nonlocal) terms, we obtain H — X]a=p a Ha + Hjy, where 



Here the sum consists of charge p and spin a chan- 
nels. Hh (X J dy cos(\/8$o-(y)) describes the undiago- 
nalizable backward scattering term^. Ha and <I>a are the 
bosonic operators satisfying the commutation relation: 
[$a(y),nQ/(2/')] = iSa,a'S{y - y') ■ The renormalized ve- 
locity and Luttinger exponents are 
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where 5e„/</.„ = i'^94.a T (2.g2,a - gi,||)) and g,^p/„ = 
±5i,-L)- Here g^^y^ = J ^ [Vi/o{qx,0)], 
92,\\/± = J ^Vi/o{qx,0)cos{2qa;kFlo), and gi,i\/i_ = 
I ^^^'^i/oilx, ^kp) are defined as the usual g-ology inter- 
action in the Luttinger liquid theory^ with kp being the 
Fermi momentum. Vi{q±) and Vb(q±) are the intra- wire 
and inter-wire interaction matrix elements, respectively. 
To simplify calculations we model the screened Coulomb 
interaction by using y(q) = (47re2A2/eo)e-lil'^o where 
eo is the static dielectric constant and Aq is screening 
length. The qualitative results obtained below should 
not be sensitive to the details of the screening potential. 

The ferromagnetic transition occurs as the spin stiff- 
ness, ujvo- = Ucr/Kcr = vp{l + (jcf,^), bccomes zeroiS, or 
(7i|| = 2TTVp+2(g4,^+g2^a-)- In general the low energy Lut- 
tinger liquid parameters should be renormalized by the 
backward scattering. Hi,, and therefore the phase bound- 
ary obtained from the bare Luttinger parameters should 
be modified also. However, when in PSFM phase, the 
spin stiffness is negative so that higher order derivatives, 
like dy^a, has to be included to stablize the system and 
to give a nonzero spin density, ps oc dy^c/^- As a result, 
the sine-Gordon backward scattering will oscillate in real 
space and hence become negligible after averaging in the 
thermodynamical limit. Therefore for simplicity we may 
assume that the renormalization effects are not very se- 
rious so that the phase boundary of the PSFM state can 
still be estimated roughly by using the bare Luttinger pa- 
rameters as stated above. The critical behavior of similar 
transition has been also discussed very recently^". 

In Fig|3we show the calculated critical inter- wire dis- 
tance as a function of magnetic field for various single 
wire electron densities, rtg. PSFM occurs in the large field 
and small distance regime. At zero distance, g2/i,a — 0, 
and therefore the critical field (wccr) is the minimum 
field strength for the backward interaction (ffiji) to be 
dominant. On the other hand, in the extremely large 
field regime, the Fermi velocity approaches zero. The 
critical distance (dcr) is now determined by the competi- 
tion between the backward scattering and the forward 
scattering in the spin channel. In large density limit 
(kpXo = mieXo > 1) we can obtain the analytic ex- 
pression of LUc.cr and dcr- ^c.cr ~ ujq\J 2r^^ e^'^^''^^'^'^ — 1 
and dcT ~ Aoe"^^'^^'^'')^ , where = rae^jtQTxkp is the 
ratio of the average potential and kinetic energies. We 
also checked explicitly that in the parameter regime we 
consider here the pseudospin polarized state is always 
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energetically unfavorable compared to the (easy-plane) 
pseudospin ferromagnetic phase. 

We now discuss how such PSFM phase can be observed 
in realistic experiments. In this phase the system has 
quasi long-range order characterized by the presence of a 
Goldstone mode. Tunneling spectroscopy has been used 
to observe similar modes in the QH bilayersS and can be 
also applied to the present system. We expect a strong 
enhancement of the tunneling conductance at small volt- 
age bias when the system enters the PSFM state. An- 
other approach to demonstrating the ID PSFM in the 
double wire system is to perform the Coulomb drag ex- 
periments. Such experiments have been done on 2DSi and 
ID^^ semiconductor heterostructures in recent years, and 
the drag resistance, i?^, is a direct measure of the effects 
due to inter-wire interactiorJ^. In the literature without 
magnetic field or interwire coherence, the drag resistance 
behaves differently in the two different regimes: In the 
perturbative regime Rd vanishes in low temperature limit 
{Rd oc <C e'^ /h)i^^; in the strong interaction regime, 
however, the backward scattering between the two wires 
becomes relevant ^■^ and opens a gap A in the energy spec- 
trum, corresponding to the formation of a locked charge 
density wave phase (LCDW) with a divergent drag resis- 
tivity Rd oc exp (A/T) in low temperature regime. 

To analyze the drag resistance in the presence of inter- 
wire coherence, it is useful to employ the Hartree-Fock 
(HF) approximation. This approach neglects long wave- 
length fluctuations present in ID systems, but we ex- 
pect these fluctuations give rise only to small correc- 
tions in the drag resistance deep inside the PSFM phase. 
The HF Hamiltonian then can be easily diagonalized by 
transforming the electron operators into the symmetric 
(cj + c| J.) and the antisymmetric (cj ^ — c| j.) channels 
with the eigenenergies, i?^ — k'^/2m* -I- =F A^ — Wq. 
Here Sfe and A^ are the intra- wire self-energy and the 
IWC gap respectively, and Wq is the shift of the band 
energy in response to the reconstruction of the ground 
state due to coupling to leads, see Fig. IHb)-(c). For 
simplicity, in our calculation we neglect the momentum 
dependenace of and A^ and approximate them by 
their values at k = 0. Within this approximation, we 
obtain (at zero temperature): 



Ao 



y^g-(rf/2AQ) 
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(6) 



where Vi = e^Xo/e^Xx, X^ 



Zq/4, and Xy 



Ag + ricoh = 2ne is the total electron density 

of both wires in the coherent regime. In above equations, 
we have assumed that all electrons fall into symmetric 
band. This is justified because the bottom of the anti- 
symmetric band can be shown to be above the chemical 
potential by /S.E = 2Ao--4?7^£'i? > 0, when the magnetic 
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FIG. 3; (a) Typical setup for conductance experiment of the 
double wire system, where the two wires interact in the middle 
regime {Q < y < L) and are connected to ideal ID reservoir 
in the left (y < 0) and right {y > L) hand sides. The upper 
(active) wire is biased by a voltage V, while the lower (passive) 
wire is biased by Vr and Vl with currents /|/| in the two wires 
respectively, (b) and (c) are the band energy for electrons 
in the incoherent reservoirs and in the coherent double wire 
regime respectively. The upper and lower bands in (b) are for 
the antisymmetric and symmetric bands respectively. 



field is large enough (Fermi energy Ep — cx B ). 

To calculate the drag resistance in a typical experi- 
mental setup, Fig. Ela), we first note that the drag 
resistance [Rd = [Vr — Vl)/I] for — 0] can be ex- 
pressed through the conductance of symmetric currents 
[G+ = /t /V for Vl = V,Vr^O and hence = /x] and 
the conductance of antisymmetric currents [G_ = /| /V 
for V/, = 0, Vr = y and hence /| = according 
to: Rd = GZ — . The symmetric and antisymmetric 
conductances, G±, in the presence of inter- wire coherence 
at temperature T can be easily derived to hei^, 
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where tg/a and r^/^ are the transition and reflection co- 
efficients for the symmetric/antisymmetric channels re- 
spectively. For simplicity we assume that Aq is constant 
for < y < L and vanishes outside this interval (the 
shaded area of Fig. |Sl[a)). We then obtain 
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{P - K^)sin(KsL), 
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where D ~ (fc^ + n'j.) sin(KsL) -f- 2ikKs cos(ksL) and 
Ks = -^/fc^ + (4?7^ — l)kp. The momentum k is re- 
lated to energy E in Eq. ^ hy E = fc^/2m*. Va 
is also given by Eq. (jHJ, replacing Kg ina, where 
Ka = v'2C- (4r?2 - l)k'j, - fc2 and C = Ao/Ep- 

At zero temperature the conductance and hence the 
drag resistance exhibit periodic dependence on the num- 
ber of electrons. At intermediate temperatures, vp/L <C 
T Ep, these oscillations are smeared out yielding 
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FIG. 4: Drag conductance as a function of ^ = Aq/Ef, follow- 
ing Eq. 0. Results for two electron densities, rj, are shown 
together. Inset: Drag conductance as a function of magnetic 
field for d = O.OSAq. rirea = 0.6 and 0.7 x 10^ cm"^ for the 
lower and upper curves respectively. 



where Rq = 27r/e^. In Fig.Qlwc show the calculated drag 
resistance as a function of <^ = Aq/ Ep. It is negative 
when ^ is small, but becomes positive with increasing ^ 
and eventually saturates at (1 — l/2r])RQ. 

When applying above results to realistic system, one 
should remember that due to the repulsive inter-wire 
interaction, the total electron density in the coherent 
regime, ricoh, should be smaller than the total electron 
density in the incoherent wires, 2nres. Such electron de- 
pletion is negligible in bulk materials due to long range 
Coulomb interaction and formation of dipole layers on 
junction surfaces. The latter ensure that bringing two 
bulk 3D materials in contact and equilibrating their elec- 
trochemical potentials does not change their densities. 
In ID systems, however, the dipole layer effects are 
greatly reduced so that the ratio of electron density in- 
side the IWC regime to the density in the reservoirs, 
T] = ncoh/2?T.i-es, may be appreciable smaller than one. 
Within HF approximation, we obtain for small di^ 
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where kp ~ Trnics is determined by the electron density in 
the reservoir. Using the same parameters as in Fig. [3 we 



plot the drag resistance as a function of magnetic field at 
a given inter- wire distance and electron density rij-cs. We 
note that a finite drag resistance {Rd does not scale with 
the wire length at T = 0) is a signature of the coherent 
state. The origin of this effect is the indistinguishibility of 
electrons flowing in the active and passive wires ((cjcj^) ^ 
0). Similar phenomenon has already been observed in the 
2D QH bilayer systems^. 

As mentioned above, without the magnetic field and 
inter-wire coherence, the ground state of the double 
wire system is predicted to be a LCDW for long-range 
Coulomb interaction with an infinite drag resistance at 
zero temperature. The effect of forward scattering could 
also be relevanfei^ at elevated (but still small compared 
to the Fermi energy) temperatures. Rd calculated in this 
scenario always increases as the inter-wire distance de- 
creases, due to the enhancement of inter-wire interaction. 
However, as we have shown in this Letter, when a strong 
magnetic field is applied, a finite Rd that does not scale 
with the wire length is expected to be observed when en- 
tering the PSFM phase. Combination of the above two 
results leads to the following overall description of the 
drag reistance: when the inter-wire distance is decreased 
from a large value (or the magnetic field is increased from 
zero) the low temperature drag resistance should first 
increase and reach a maximum value around the phase 
boundary (Fig. ^ and then begin to decrease to almost 
zero due to IWC when entering the PSFM phase. Such 
nontrivial behavior of drag resistance could indicate a for- 
mation of ID pseudospin ferromagnetism in small inter- 
wire distance or large magnetic fields. 

To summarize, we have shown that in the presence of 
a strong magnetic field the electronic system can become 
(pseudospin) ferromagnetic in the double quantum wire 
system. We further demonstrat that the low tempera- 
ture drag resistance has a non-monotonic behavior near 
the phase transition boundary, which should become ob- 
servable in the present or near future experiments. 

We appreciate fruitful discussion with S. Das Sarma, J. 
Eisenstein, B. Halperin, H.-H. Lin, Y. Oreg, M. Pustilnik, 
A. Shytov, A. Stern, A. Yacoby, and M.-F. Yang. This 
work was supported by Harvard NSEC and by the NSF 
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